Orbital Zeeman effect: Signature of a massive spin wave mode in ferromagnetism 
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By deriving the quantum hydrodynamic equations for an isotropic single-band ferromagnet in an 
arbitrary magnetic field, we find that a massive mode recently predicted splits under the action of 
the field. The splitting is a peculiarity of charged fermions and is linear in the field to leading order 
in q bearing resemblance to the Zeeman effect in this limit, and providing a clear signature for the 
experimental observation of this mode. 
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Magnetism in solids has been one of the most studied 
subjects in physics for the past decades. In particular the 
behavior of propagating spin waves in metallic materials 
has now a long list of contributions in the literature. Spin 
waves in ferromagnetic materials were first predicted by 
Bloch J!]] and Slater 0j and later observed in iron by 
Lowde. || These first theories were constructed based 
on lattice models of local moments. Predictions for spin 
waves using Fermi liquid theory [0 were first made by 
Silin for paramagnetic systems. (5T Paramagnetic spin- 
waves were found to propagate only under an applied 
magnetic field. Abrikosov and Dzyaloshinskii || were 
the first to develop a theory of itinerant ferromagnetism 
based on Landau's theory of Fermi liquids. Although 
correct at the phenomenological level, the microscopic 
foundations for this theory were estabilished only much 
later by Dzyaloshinskii and Kondratenko. Q| 

In this paper we derive the hydrodynamic equations of 
the ferromagnetic Fermi liquid theory (FFLT) for a finite 
magnetic field and show that an intrisic degeneracy of a 
recently predicted massive mode || exists and is lifted 
under the external field, similarly to the Zeeman split- 
ting of a single spin in a magnetic field. Underneath this 
effect we identify the breaking of chiral symmetry by the 
external magnetic field in the case of charged fermions, 
which relates to the orbital motion developed under the 
Lorentz force. Besides being a clear signature of the mas- 
sive mode, it is suggested that quantization of this mode 
will lead to a system of massive magnons with "up" and 
"down" states that split under an external field. We es- 
timate the values for the fields where the effect will be 
observed in typical weak ferromagnets. 

For isotropic metals, the equations describing dynam- 
ics of paramagnetic spin waves are similar to the ones 
that result from FFLT in the small moments limit. This 
is expected since in this limit quasi-particles can be de- 
fined and one recovers the kinetic behavior of ordinary 
Fermi liquid theory. However, FFLT rests on the assump- 
tion of a quite different, symmetry- broken, ground state 
and the resulting spin waves propagate with no external 



magnetic field present. The Goldstone mode associated 
with spontaneously broken spin rotation invariance has 
been the paradigm of spin-waves in an isotropic single 
band ferromagnet and can be derived from FFLT (as 
well as from lattice models). However, FFLT contains 
spin-wave modes that have not yet been observed, as 
pointed out in Ref . ||] , where the proper hydrodynamics 
and parameters for the propagation of the lowest in en- 
ergy of these modes have been studied. Such mode is not 
a Goldstone mode, hence its dispersion may in principle 
be gapped for q = 0. Indeed, for low q and zero external 
magnetic field, its dispersion has been shown to be of the 
form u> = uj^ — aq 2 , j^] and propagation is possible only 
in the quantum hydrodynamic regime (at temperatures 
for which collisions are almost absent). This is in con- 
trast to the Goldstone mode which is purelly quadratic 
and also propagates at collision dominated temperatures. 
JTo| Propagation of the massive mode also requires that 
the interactions have some finite amplitude with p-wave 
symmetry. Lack of observation of such a mode in neutron 
scattering experiments can be attributed to the low spec- 
tral weight the mode carries at small q (~ q 2 when rated 
against the Goldstone mode), however indirect signs of 
its existence have already been observed. |sj 

We begin our derivation by writing the equal 
time propagator for a system of interacting fermions, 
G ai( y 2 (i"ir 2 ) = (aj Ti (ri)a (T2 (r 2 )), where a is the usual 
fcrmion anihilation Heisenberg operator and the aver- 
age is taken with a ferromagnetic ground state. The 
equation of motion for G is, in spin space, c*tG(rir 2 ) + 
(i/n)/dr'[W(nr')G(r'r 2 ) - G(nr')H(r'r 2 )] = g[G], 
where TL is the interaction Hamiltonian and the func- 
tional g[G] is the full time-derivative of G. Following an 
usual script, one turns to a mixed r - p representation by 
defining « p (r) = / dr' e~ lp r ' / h H(r + r'/2, r - r'/2), and 
likewise for A/" p (r), which brings the equation of motion 
to the form 

oW P (r)+ J Dt [?V +rik/2 (r' + ?Ls/2)AV(r') 

-N- p ,(r')H p ^ hk/2 (r' -hs/2)] = J[Af p ], (1) 
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where Dt = idr'dsdp'dke l tt p '- p > s+( - r - r ') ^/{2ir) 6 h. Ex- 
panding Eq.([j]) in s and k yields a formal series in H, 
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where we have (anti) commutation if m + n is (odd) even 
and sumation over repeated Greek indices. 

It is clear from this expression that the repeated ac- 
tion of the derivatives on the fluctuating fields generate 
coefficients that are higher order in q for higher powers 
of H. If we keep terms up to first order in q, 

d t M p (r) + % - [H p (r) ,A/p (r)] - \ {W p (r), A/p(r)} 

+ ^{^ P (r),Wp(r)} = J[^ p ], (3) 

where square braces indicate the usual commutator while 
curly braces are Poisson brackets. Therefore, keeping 
first order in q leads to the same result obtained by drop- 
ping terms in the series that are explicitly 0(h 2 ). How- 
ever, H p (r) is itself a function of U, which disables, at 
this point, comparison with a semi-classical approach. 

For the ferromagnetic Fermi liquid there are two dis- 
tinct underlying scenarios that we call (following Ref . || ) 
classical and quantum spin hydrodynamics. We want 
here to discuss in a general way the origin of these two 
regimes. For this purpose, it should be recalled that the 
equilibrium density matrix for a ferromagnetic state is 
written in the quite general form, j(| 
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where t are Pauli matrices and mo is the equilibrium 
magnetization density, which is a conserved quantity of 
the ferromagnetic Hamiltonian. From this it is clear that 
the commutator in ([|) is zero unless there are fluctuations 
about equilibrium. 

At high temperatures, thermal fluctuations dominate 
and this commutator remains irrelevant. The remaining 
"classical" terms lead to the known Bloch-like hydrody- 
namics for the Goldstone mode. This regime is refered 
to as "classical spin hydrodynamics." 

At low temperatures quantum fluctuations dominate, 
and the commutator in ([|) becomes important. It re- 
mains finite as long as the fluctuations contain, besides 
the usual Goldstone-mode, contributions from non con- 
served quantities as it is the case of spin-current, whose 
oscillations give rise to the gapped mode studied in Ref. 
H , where this regime has been called "quantum spin hy- 
drodynamics." 

The crossover temperature for these two regimes 
has been found to scale with {uj^) 1 / 2 , where lj^ = 



2|m |(F 1 a /3 - F$)/hN{Q) is the gap, F^ a)s are the usual 
dimensionless spin (anti) symmetric interaction parame- 
ters, and N(0) is the density of states at the Fermi sur- 
face. 

The preceding discussion is general and establishes 
the origin of the massive mode. We now turn to the 
small moments limit, where kinetic equations can be 
amenably derived. In this limit, we write the effec- 
tive quasiparticle Hamiltonian and density matrix re- 
spectively as H P aa'(r) = e P (r)(5 CTCT / + h p (r) • r CTff /, and 
A/pcrcr' (r) = n p (r)S aa i + er p (r) • r aa i . These forms for the 
Hamiltonian and density immediately give 



[Hp(r),A/p(r)] = -(<x p (r)xh p (r)) 



(5) 



Here the internal field is given by 
Mr) 



?H + ^£rW<P-P'K>(r), (6) 

p't 



where Pg are Legendre polynomials, H is an external 
magnetic field, and all coupling constants are h inde- 
pendent. We see then from Eqs.(||) and (||) that the 
commutator of Eq.(|^) gives two contributions: a zero-th 
order term in Ti which is just the (Larmor) precession 
of the internal magnetization and a first order term in 
h that is the precession about the field generated inter- 
nally. Carefull examination of Eq.(|J) shows that this 
term is the only first order term in fi of the series. That 
is, provided that the effective Hamiltonian is of the Fermi 
liquid type, the long wavelength limit is equivalent to a 
Wentzel-Kramers-Brillouin or (eikonal) semi-classical ap- 
proximation. 

We want now to have the Lorentz-force operator ex- 
plicitly factorized. This is achieved by rewriting Eq.(|^) 
in a gauge transformed "frame" for which A = 0. Also, 
to leading order in the fluctuations, the non-commuting 
parts of the Poisson brackets can be dropped, giving 



V,W + s [Hp(r),jV p (r)] + {AT p (r),W p (r)) 



dp 



dp 



We see that the gauge boost factorizes a non-chiral term 
whose amplitude is proportional to e|H|. This is the con- 
tribution from the orbital motion of charged quasiparti- 
cles, a result similar to the one found in the theory of 
normal metals, [Tl| ] where it is known simply to shift the 
paramagnetic resonances. In the isotropic ferromagnetic 
case we are studying here this term plays a more essen- 
tial role. Equation (Q) is formally identical to the one 
obtained in a normal metal. The difference rests on the 
broken symmetry ground state, whose density is given by 
Eq.(@). 

Let TZ + ( - 1 be a simultaneous proper (improper) rota- 
tion of real and momentum axes and likewise SU + ^ for 
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the spin axes. Examination of Eq.(|7|) reveals that in the 
absence of a magnetic field, IZ^ are symmetry operations 
both in the paramagnetic and ferromagnetic cases, SIA~ 
is not a symmetry in either case, and SIA + distinguishes 
the two cases, being a broken symmetry in the ferromag- 
netic case. The presence of H adds two terms to the 
equation, one that couples to spin through h p and the 
orbital one discussed above which couples to charge. In a 
charged paramagnet, the former will break SIA + symme- 
try yielding propagating paramagnons as a consequence, 
while the latter will break 1Z~ (chiral) symmetry. How- 
ever, since the spin modes are sustained by a finite H, 
the degeneracies associated with chiral symmetry are 
lifted with H, before spin waves can propagate. 

In the isotropic ferromagnet, spin waves can propa- 
gate in the absence of H, for the SU + is a spontaneously 
broken symmetry and the internal "sustaining" field is 
provided by m . Hence propagating modes exist before 
the formation of orbits with the breaking of 1Z~ by a 
finite H, and some of these modes may be degenerate. 
Breaking chiral symmetry should then lift such degen- 
eracies. Note that the Goldstone mode should not be 
affected by H (apart, of course, from shifting) for the 
term breaking chiral symmetry only couples to charge. 
However, the gapped spin-wave mode splits as it is seen 
in the dispersion relations of Fig. ([!]). 

In order to see how this happens, we finish linearizing 
Eq.(^), and then solve it in the hydrodynamic limit. ||] 
In particular, since we are seeking equations on the total 
magnetization density m = tr[T^ p 7V p ], and spin cur- 
rent tensor j rri = tr["r^ p (cW p /<9pi)A/p], we trace the 
product of r with Eq.(^|) and keep only terms that are 
linear in the fluctuations. The result is a set of coupled 
equations, 

d t 5m + V a j tTa = -7H0 x 8m - 7<5H x m , (8) 

which is just the continuity equation, and 

dtja a - c 2 s \7 a (Sm + xo^H) + T3*e az pj„p 

= -7(H -w+m ) x ja a -^* D ja a - (9) 

Here Sm = m — mo, (5H = H Ho where Ho is a uniform 
magnetic field applied parallel to the z axis defined by the 
equilibrium magnetization mo, xo = 7?jAT(0)/2|l + Fq\, 
H* = e|H |(l + Ff/3)/m* is the orbital amplitude, c 2 s = 
0173) 1 1+^1(1+^/3) is the squared spin wave velocity, 
and u* D = (1 + Fi/3)/td is the p-wave amplitude of the 
scattering integral for a spin-diffusion relaxation time td- 
We consider the response to a driving field SH which is 
transverse to z and oscillates with frequency ui. It is then 
easy to see that the three longitudinal components z-j aa 
vanish together with 8m z . This means that there are 
no magnetization gradients in the longitudinal direction 
and there is no transport of longitudinal magnetization 
in any direction. The remaining 6 components can be 
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FIG. 1. Dispersion relations at vanishing temperatures. 
H* = e|H |(l + Fo73)/m*. 

conveniently "folded" onto the complex "vector" j + = 



w ■ Jo 



which measures the transport of 



idrriy under the 



transverse magnetization 8m + = 5m 
driving field 5H+ = 5H x + iSH y . Solving Eqs~.(§) and (|) 
for the Fourier components yields [fL2[ 



SH+c 2 s 



and 5m 



[T I m o I -Xo(w-wl)] 

x [f2 2 q — iSliL* x q - (H* • q)H*] , (10) 

-5H + -/\m \P 3 (Q)/Pi(n), where fl = 
lo — ljl — zq, ojl is the Larmor frequency, zq = 

lo+ - ilu* d , p 3 (fi) = n 3 + (c 2 s q 2 Xo /i\mo\)n 2 - w 2 n - 

( C 2 Xo /7|m |)(H* ■ q) 2 , P 4 (0) = tt 4 + z Q n 3 + (c 2 q 2 - 
R* 2 )n 2 - z H* 2 n - c 2 (H* • q) 2 , and H* = Wz. For 
simplicity we specialize to the case Ho _L q (indicated 
by 9 = 90° in the figures). This condition is obtained in 
practice by sheding neutrons or photons adequately on 
a finite sample's surface, and is the condition that max- 
imizes the non-chiral orbital contribution in Eq.([l0|); in 
the infinite system, we can think of transverse magnons 
if we wish to quantize the hydrodynamics. Longitudinal 
magnons will be treated elsewhere. Jl2| In Fig. (Q) we see 
the dispersion relations for different values of the external 
field. These are obtained by looking into the free modes 
(<5H = 0). It is clear that the degeneracy is only exact 
at q = 0. We see also that it constitutes a three-fold 
degeneracy, however, the two constant branches shown 
(uj — lol = Acj = uj^) are spurious when Ho = in 
the sense that there is no spectral weight associated with 
them, as seen in Fig.(||) where S(q,u>) tx Im[5m + /5H + ] 
(which is the important quantity in experiments) is plot- 
ted as a function of uj and H*. Under a finite magnetic 
field only one of the spurious branches develops weight 
while the degeneracy is lifted. 
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FIG. 2. Spectral weight in arbitrary units. The splitting 
is clearly seen as the magnetic field increases. Numbers on 
the upper left corner give the amplitude of one plot relative 
to the next. These amplitudes are of order of 10 -2 the (not 
shown) Goldstone mode amplitudes which do not split under 
the external field. The curves are translated vertically to fa- 
cilitate analysis. The values of H* correspond to hundreds of 
Gauss in, e.g., MnSi. 



the same way as in the Zeeman effect of a single spin. The 
possible existence of (low temperature) magnons with 
mass and a Zeeman-like degeneracy makes room for new 
questions realated, e.g., to the collective behavior of a 
handful of these excitations. Before these questions are 
put forward it is, however, prudent to focus on the issue 
of whether these results can be observed by traditional 
experiments. The spectral weight of this mode, of or- 
der of 10~ 2 the usual spin waves (see Fig.(||)) stands as 
an obstacle. Our main emphasis here is whence on the 
additional criterium such a degeneracy provides to track 
down the massive mode. In a companion longer article 
we also discuss how these results are expected to show in 
a conduction ressonance experiment. |fl2|| 
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The important point to note is that the frequencies 
of the Goldstone mode (the traditional spin- waves) only 
shift their loci as the external magnetic field increases 
whereas the gapped mode splits. This can be seen in 
Fig. (||) where the behavior of these dispersions with the 
external field is shown (including the Goldstone mode). 
It is then clear that the splitting of the gapped mode 
should be considered as a distinctive feature in experi- 
ments searching for a direct observation of it. It should 
be helpful to make a statement about the quantities in- 
volved for some traditional materials: From data found 
in the literature we estimate the fields shown to be in the 
range of hundreds of Gauss for typical weak ferromagnets 
like MnSi, ZrZn2, and M3AI. This corresponds to a gap 
of tenths of meV in these materials. 1131 
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FIG. 3. Linear behavior of the dispersions with H*. The 
effect of H* on the Goldstone mode (ui — ujl ~ 0) is seen to 
merely shift its frequency. 

For small values of H* and q (up to O.lwj 1 "), the dis- 

„2„2 /, ,+ 



persions may be put in the simple form Aw 



<z7< 



for the Goldstone mode and Alu = ±H* + u>^ — aq for 
the gapped mode. |l4|] Quantization of the theory in this 
limit will yield massive magnons with an "up-down" de- 
generacy which is lifted by the magnetic field, in much 
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